Sigma term and strangeness content of octet baryons 
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Abstract 

By using lattice QCD computations we determine the sigma terms and strangeness content 
of all octet baryons by means of an application of the Hellmann-Feynman theorem. In addition 

■ to polynomial and rational expressions for the quark-mass dependence of octet members, we 
{SJ ' use SU(3) covariant baryon chiral perturbation theory to perform the extrapolation to the 

■ physical up and down quark masses. Our Nf = 2 + 1 lattice ensembles include pion masses 
' down to about 190 MeV in large volumes (M 7r L>'i), and three values of the lattice spacing. 

P-H ■ Our main results are the nucleon sigma term g^n = 39(4)(^ 8 ) and the strangeness content 

(T) . un = 0.20(7)(1}7). Under the assumption of validity of covariant baryon xPT in our range 

of masses one finds = 0.276(77)(lg2)- 

'+->'r 

PJ 1 Introduction 

<D 

The nucleon sigma term, and the nucleon strangeness content are phenomenologically im- 
fS| 1 portant quantities that are not directly accessible to experiment. They are related to tt — N 

and K — N scattering lengths, to the quark-mass ratio m s /m U( i (where m M d = (m u +m c [)/2), 
to the hadron spectrum and even to counting rates in Higgs boson searches. They also play 
" a key role in the direct detection of dark matter (DM): sigma terms relate the effective DM - 

nucleon coupling to the fundamental Lagrangian parameters that couple the DM particle to 
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ON I the quarks. A precise first principle determination of these quantities is thus very interesting. 

The sigma terms are defined as 



X 



<JnN = (N(p)\m ud (uu + dd)(0)\N(p)) (la) 
a- ssN = 2(N(p)\m s ss(0)\N{p)). (lb) 

The usual way of computing a n N is by using tt — N scattering data. a T N cannot be directly 
obtained in this way, but available data can be extrapolated to the Cheng-Dashen point, 
although the machinery involved in this determination is complicated. xPT can then be 
used to reliably extrapolate from the Cheng-Dashen point to the chiral limit [HE]. 

Such an analysis was carried out in the 1980's [HE], resulting in a value of o-kN = 45 ± 8 
MeV. A later study claims to obtain a higher value [3] , a^N = 64 ± 7 MeV which is almost 
2 standard deviations away. Also it has been pointed that the uncertainties associated with 
these determinations are already affecting the interpretation of (direct) dark matter search 
experiments [5] , although lattice simulations are already helping to change this situation [5] . 

Sigma-terms of other octet members are also of phenomenological interest. Recently they 
have been used in the context of the hadron resonance gas model for estimating quark-mass 
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effects in QCD thermodynamics calculations [6j. Moreover, as we will see, an analysis of all 
octet members can be further used to constrain the strangeness content of the nucleon. 

In this paper we compute the sigma terms of the octet baryons a- w x-,&ssX for X = N, A, S, 5. 
The Feynman-Hellman theorem applied to QFT relates the sigma terms to the dependence 
of baryon masses with respect to quark masses through: 

dM x 

o-kX = m ud~R (2a) 

(Jssx = 2m s -^—. (2b) 

This opens the possibility of computing the octet sigma terms via the variation of octet 
masses with respect to quark masses in a lattice simulation, which is the approach that we 
use here. A preliminary account of this work has been given in [7]. There are other interesting 
quantities directly related to the previously mentioned sigma terms that are also useful for 
phenomenology 

2(X(p)\ss(0)\X(p)) 
yX (X(p)\(uu + dd)(0)\X(p)) [ > 

f m ud (X(p)\(uu + dd)(0)\X(p)) 

fudx = Wx (3b) 

m s (X(p)\ss(0)\X(p)) 

fsX = w x • (3c) 

The first is called the strangeness content, and the other two are sometimes referred to as 
" dimensionless sigma terms". The latter two are directly related to the conversion of the 
fundamental DM-quark coupling to the effective DM-nucleon coupling. 



2 Simulation details and ensembles 

The gauge and fermionic actions, as well as the algorithms used are described in [8l[9]. Here 
it suffices to mention that we simulate QCD with two degenerate light quarks and a heavier 
strange quark, and that we use tree level improved Wilson fermions. 

To set the lattice spacing corresponding to each value of the coupling (3 and fix the quark 
masses to their physical values, we use M n , Mk, and M^. We extrapolate the values of 
aM n ,aMK,aMQ to the point where the ratios aM^/aM^ and clMk /clMq agree with the 
experimental values^. 

As shown in Table [I] at /3 = 3.3 and (3 = 3.7 the strange quark mass is held fixed, whereas 
for (5 = 3.57 we simulate at three different values of m s to have some lever arm to perform 
the small extrapolation to mf ys (see also Fig. [1]). Our data sets cover a wide range of pion 
masses from M n ~ 190 MeV up to M n ~ 680 MeV, although in this analysis we only use 
ensembles with M n < 550 MeV. 

On every ensemble we measure the octet masses Mn,M\,M-£ and M= with valence quark 
masses equal to sea quark masses (only the unitary theory is considered). It is worth men- 
tioning that these ensembles have previously been used to accurately predict the light hadron 
spectrum [9], including the masses of the octet baryons. 

2 Experimental inputs are corrected for isospin breaking and electromagnetic effects according to [S]. 
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arriud 


am s L 3 xT 


Trajectories 


aM n 


aM K 




-0.0960 - 


0.057 16 3 x32 


10000 


0.4115(6) 


0.4749(6) 




-0.1100 - 


0.057 16 3 x32 


1450 


0.322(1) 


0.422(1) 


3.3 


-0.1200 - 


0.057 16 3 x64 


4500 


0.2448(9) 


0.3826(6) 


-0.1233 - 


0.057 24 3 x64 


2000 


0.2105(8) 


0.3668(6) 




-0.1233 - 


0.057 32 3 x64 


1300 


0.211(1) 


0.3663(8) 




-0.1265 - 


0.057 24 3 x64 


2100 


0.169(1) 


0.3500(7) 




-0.0318 


,-0.010 24 3 x64 


1650,1650 


0.2214(7),0.2178(5) 0.2883(7),0.2657(5) 


3.57 


-0.0380 0,-0.010 24 3 x64 


1350,1550 


0.1837(7),0.1778(7) 0.2720(6),0.2469(6) 


-0.0440 


,-0.007 32 3 x 64 


1000,1000 


0.1348(7),0.1320(7) 0.2531(6),0.2362(7) 




-0.0483 


,-0.007 48 3 x 64 


500,1000 


0.0865(8),0.0811(5) 0.2401(8),0.2210(5) 




-0.007 


0.0 32 3 x96 


1100 


0.2130(4) 


0.2275(4) 




-0.013 


0.0 32 3 x96 


1450 


0.1830(4) 


0.2123(3) 


3.7 


-0.020 


0.0 32 3 x96 


2050 


0.1399(3) 


0.1920(3) 




-0.022 


0.0 32 3 x 96 


1350 


0.1273(5) 


0.1882(4) 




-0.025 


0.0 40 3 x96 


1450 


0.1021(4) 


0.1788(4) 



Table 1: Parameters of our simulations. The errors quoted here are purely statistical. These results 
correspond to one of the 18 two-point function, time fit intervals that we use in our estimate of systematic 
uncertainties. In this particular analysis, the scales at /? = 3.3, 3.57, 3.7 are a -1 = 1616(20) MeV, 2425(27) 
MeV, 3142(37) MeV, respectively. 



3 Chiral extrapolation 



As already stated, we simulate QCD for values of the light quark masses larger than the 
physical values, whereas the strange quark is close to its physical value with some lever arm 
to perform a small extrapolation. Thus our data requires an extrapolation in the light quark 
masses and an interpolation in m s . 

For this, we need to describe the quark-mass dependence of octet baryons. We will always 
use the tree level SU(3) chiral relation to express the quark-mass dependence through the 
meson mass dependence via 

m ud oc Ml 

m s cx Ml = 2M 2 K -Ml. (4) 



To study this mass dependence and the associated model uncertainty, we consider two very 
different approaches [9l|10]. First we consider a regular expansion in quark masses around a 
non singular point where none of the quark masses vanish. For this purpose we define the 
expansion variables 



A w 



Ml 


- (M£ en ) 2 




Ml 


Ml 


- (M£ hys ) 2 


Ml 



(5) 



where (Mf n ) z = \ 
included in the fit. 



M? hys Y 



+ {M^y 



and M£ denotes the mass of the heaviest pion 
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Figure 1: Overview of our simulation points in terms of M n and y^2M^ — M%. The former gives a 
measure of the isospin averaged up and down quark mass while the latter determines the strange quark 
mass. The symbols refer to the three lattice spacings, and the physical point is marked with a cross. Error 
bars are statistical only. 

Note that changing the expansion point in polynomial-like formulae is simply a reshuffling 
of the coefficients. When interested in checking the SU{3) symmetric line m u d = m s it is far 
more convenient to define our expansion variables as 



where the SU(3) symmetric line is characterized by A^ = A Ss . Because the extrapolation in 
m s is small, a linear term in Ag S suffices 

Second we consider a singular expansion around the SU(3) chiral point m U( i = m s = 0, which 
allows for a much more constrained expansion. Low-energy processes in QCD have been 
extensively studied in the framework of chiral perturbation theory. The seminal work |11] 
and its success in explaining meson observables, have pushed the study of baryons within 
the same framework. However baryon %PT is far more involved. The main difficulty is 
that baryon masses Mb are not small in comparison with the scale of chiral symmetry 
breaking (A x ~ 1 GeV), and Weinberg's power counting theorem [12] fails: higher-order 
loop corrections contribute with powers of Mb /A x and are no longer small. This results in a 
very slow (practical) convergence of the series [13] . This convergence is partially improved by 
treating the baryons as heavy degrees of freedom, in what is known as the heavy baryon chiral 
perturbation theory [II]. Although it is reasonable for dealing with the mass dependence 
of baryons with light quarks, experience seems to show that observables that depend on the 
strange quark are not well described in this framework (see for example |15H17j ). A number 
of possibilities have been proposed to improve the convergence and include the strange quark 
in the analysis. The various flavors of cutoff B^PT (see [18] and references therein), or 
covariant B%PT |19j are some examples. In this work we explore this last possibility. The 
interested reader can consult the review [20] and references therein. 




(6) 
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3.1 Regular expansions 



We have several possibilities if we decide to expand around a regular point. We can treat 
any octet member X = N, A, E,S as independent and use a Taylor expansion. Using the 
variables defined in Eq. (|5|) our results are well described by the ansatz 

M x = M x + a x A^ + a^Al + /3 X A Ss . (7) 

One can also use a Pade like functional form 

M _ Ml (R \ 

MX l-ofA,-a^-^A s -; W 
In these expansions, M x ,a x 2 and ft X are the 16 fitting parameters. 

In principle all octet masses should be degenerate along the SU (3) symmetric line m u d = m s . 
We can choose to impose this constraint in our functional form. In fact our data is well fitted 
by the following SU (3) symmetric regular expansion 

M X = M + a x A n + a x A 2 n + (d - of) A s - S + (C 2 - a x ) A n A Ss , (9) 

or the corresponding Pade like ansatz 

MX = 1 - a x A w - a x Al - (d - a? ) A s - S - (C 2 - a*) A^A- SS " 

Now the fitting parameters have been reduced to 11: Mo, a x 2 and the coefficients C\ and C2 
that give the dependence of octet masses along the SU(3) symmetric line m u d = m s . 

Note that Taylor- like and Pade- like ansatze differ in higher-order terms in A n and Ag S . 
Thus the difference in physical results obtained by using these two functional forms measures 
higher-order contributions to the mass expansion. 



3.2 Covariant B^PT 

Details of the quark-mass dependence of the octet members in SU(3) baryon xPT can be 
found in Appendix[Al Here we run the index X = N, A, E, H over octet members, and index 
a = tt,K,t] over the mesons. 

Octet masses to NLO involve the nonlinear function 



/ / X \ ^ X X 

h(x) = — — j < \ 1 — I - ) arccos - + - log x > . (11) 



The masses of the octet baryons are given by 

M X = M - 4c n xM* - Ac Ss xMl + £ ^M*h (^A + M^M* + 4d Ss A^, (12) 

a=-K,K,r) a 

where c^x-.CssX are functions of the three low-energy constants (LEC), 6o)^d)^/ ; and the 
9aX can be written in terms of the axial coupling qa and the ratio of couplings £ (as shown 
in Tabs. [2] and [3]). Finally d^^dgs parametrize higher-order corrections. 
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3(l+0 2 ^ 


(1+? 2 ) „2 

(i+0 2 ^ 


(5+6£+90) 2 
6(1+^)2 </a 
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1 n 2 


1 n 2 


(1+30 2 -2 


12(l+0 2y ^ 


3(l+0 2 ^ 




12(l+0 2 ^ 



Table 2: Meson-loop couplings as a function of the ttN coupling g& and the quantity £. In the chiral 
limit gA = D + F and £ = F/D, as described in Appendix [Al 
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N 


A 


£ 




C X 


-2(2b + b D + b F ) 


-4(6o + b D /3) 


-4(6o + M 


-2(26 + 6 D - V) 


ss 
C X 


-2(bo + b D - b F ) 


-2(6o + 46 D /3) 


-26 


-2(bo + bp + 6jr) 



Table 3: Leading order octet quark-mass dependence, c n x and c Ss x, as a function of the LEC bo, bo, 6_f 
(see Appendix [A] for more details). 

To the order at which we are working only two pseudo-Goldstone boson masses are linearly 
independent, being related through the Gell-Mann-Okubo relation 3M^ = 4M^ — M%. 

When fitting lattice data to the chiral formula, it is important not to break the SU(3) 
symmetry built into this expression. In that sense it is important that all octet masses 
become degenerate when M Ss = . This rules out the possibility of fixing the meson decay 
constants F^^,^ to their physical values. Here we choose to fix all of them to the common 
value i ? | >hys = 92.2 MeV. This is correct to the order at which we are working. Note that 
we have repeated the fits including the NLO terms for the decay constants and obtained 
results compatible within statistical errors. We have also confirmed that the quality of fits 
decreases if a non SU(3) symmetric ansatz is imposed by fixing each meson decay constant 
to its physical value. 

The axial coupling gA is well known from phenomenology. The most precise value at the 
physical point is gA = 1.2695(29) [2TJ, and is expected to be close to the value in the chiral 
limit D + F. The ratio of couplings £ is not well determined experimentally, but there are 
two preferred phenomenological scenarios: £ = 2/3 and £ ~ 0.5 |22| . 

It makes sense to try our xPT fits both fixing (gA, £) = (1.2695, 2/3) or allowing them to be 
free. Regarding the higher-order contributions given by d n ,ds S , we can also choose to either 
include them in the fit, or not. In terms of how well our data are described (i.e. fit quality) 
the four possibilities are equally reasonable options. In Fig. [21 we show one possible xPT fit, 
where we have fitted gA ,£ arid we have not included the higher-order terms in the fit. 

4 Cutoff, finite volume effects and excited-state contributions 

We combine our chiral extrapolation with the continuum extrapolation. It has been demon- 
strated in different contexts [8,9,23,24] that a smeared clover action is very close to be 
nonperturbatively 0(a)-improved with cutoff effects being 0(a 2 ). Nevertheless we cannot 
rule out the possibility that our results have linear discretization effects. 

Guided by our experience [9] in determining light hadron masses we parametrize cutoff effects 
with the substitution 

M x ^ M x (l + c x a p ) (13) 
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Figure 2: Example of a xPT fit corresponding to a particular choice of fitting interval for the correlators. 
In this particular case we choose to fit gA and £, and not to include higher-order terms. We have only 
used data with < 410 MeV. The correlated \ 2 — 38.5 for 34 degrees of freedom, yields a fit quality 
of - 0.27. 



with p = 1,2. The cutoff effects are small enough in our data that they make it difficult 
to distinguish a and a 2 . In fact their coefficients are, for the most part, compatible with 
zero within statistical errors. Thus we will also include fits without any cutoff corrections 
altogether. 

Stable hadron masses have leading finite- volume effects oc e~ MnL [25J. In all of our ensem- 
bles the bound M 7r L>4 is maintained, implying negligible finite-volume effects. The detailed 
analysis of jS], with additional ensembles in smaller volumes, shows that finite- volume cor- 
rections are below the statistical accuracy of our data. Thus we do not add finite-volume 
corrections to our analysis here. 

To estimate the possible contamination by excited states in the extraction of the masses from 
the correlators we use 18 fitting time intervals: t m i n /a = 5 or 6 for (3 = 3.3, t m \ a /a = 7, 8 or 
9 for = 3.57 and t min /a = 10, 11 or 12 for /3 = 3.7 [9]. 



5 Determination of the uncertainties 



Our complete analysis includes a total of 8 formulae to extrapolate to the physical mass point: 
Four of them are regular expansions either imposing or not the SU(3) flavor constraint, and 
either in the Taylor-like form or Pade-like form. The other four functional forms are derived 
from SU(3) covariant baryon xPT, where we can choose to parametrize or not higher-order 
contributions, and to fit gA and £ or fix gA to its physical value gA = 1.2695(29) [21] and £ 
to a reasonable phenomenological valu^l £ = 2/3 [22] . 

To have more control over the uncertainties associated with higher-order terms in the mass 
extrapolation, we impose two different pion mass cuts: M n < 410 MeV and < 550 MeV. 

As explained in Sec. 0] we consider three possibilities to parametrize cutoff effects which are 
compatible with our data: either assume that they are absent, or parametrize them as 0(a) 
or 0(a 2 ). Finite- volume corrections are not considered because they are below our statistical 

3 We have checked that using £ = 0.5 (the other common phenomenological scenario) leads to very similar 
results. 
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tributions of different pion mass cuts. tributions of different functional forms. 

Figure 3: Distribution of values for ovjv in background (yellow), and how different pion mass cuts (left) 
and functional forms (right) shape it. 



accuracy. Finally we repeat the full analysis with 18 different time-fitting intervals for the 
correlators to estimate excited state contributions. 

This strategy leads to8x2x3xl8 = 864 different procedures to estimate the physical value 
of each of the sigma terms of the octet members. The data are fitted (taking correlations 
into account) using all of the previously explained procedures. Each of the 864 results is 
weighted with the fit quality |j| (p- value or Q- value) to produce a distribution of values. The 
median (typical result of our analysis) is taken as our final value. The 16— /84— percentiles 
(i.e. the values which, in a Gaussian distribution, correspond to ±lcr deviations from the 
median) yield the systematic uncertainty of our computation (see Fig. [3J). This systematic 
uncertainty measures how different parametrizations for mass extrapolation, cutoff effects and 
excited-state contributions affect the final result. The statistical uncertainty is determined 
by bootstraping the whole procedure 2000 times, and computing the variance of the medians. 

It is important to note here three points: 

• Most of the previously mentioned fits are "good." Our average fit quality is 0.34, and 
this does not vary significantly over our different procedures. For example fits with 
M n < 410 MeV have an average fit quality of 0.38 whereas for M w < 550 MeV it is 0.29 



(see Fig. 3(a)). For covariant BxPT, Taylor and Pade functional forms the p- values 



are respectively 0.30,0.47 and 0.29 (see Fig. |3(b)j ). The fit quality is almost insensitive 
to how we choose to parametrise the cutoff effects, or what time intervals we use to fit 
correlators. Although there is some variation, we cannot rule out any of our fits based 
on how well our data are described. 

Adding more variations does not increase the systematic uncertainty. For example 
adding a pion mass cut of M n < 680 MeV (i.e. including all our ensembles in the fit), 
the covariant BxPT functional form gives fit qualities of 2 x 10 -6 . Even if they are 
included in our analysis, these fits do not contribute to our distribution of values or to 
our final results and estimates of systematic uncertainties. 

On the other hand, eliminating some of the analysis typically results in a compatible 
final value but with a smaller systematic uncertainty. For example had we chosen to 



4 The confidence of fit is defined as Q = f£ dz P{z,d) where P{z,d) = z d/2 ~ 1 e' z/2 /[2 d/2 r(d/2)] is the 
probability distribution function to obtain y 2 = z in a fit with d degrees of freedom. 
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<r*x [MeV] 


(j s - sX [MeV] 


yx 


fudX 




N 


39(4)(tf) 


67(27)(t£) 


0.20(7)(ii?) 


0.042(5)(+f ) 


0.036(14)(i^) 


A 


29(3)(^) 


180(26)(+$) 


0.51(15)(±f 7 ) 


0.027(3)(i; o ) 


0.083(12) (j:^) 


£ 


23(3)(lf) 


245(29)(1^) 


0.82(21)CtS) 


0.019(3)(+f ) 


0.104(12)( + 3?) 




16(2)(+«) 


312(32)(t£) 


i-7(5)(^: y 7 ) 


0.0116(18)(t£) 


0.120(13)(t£) 



Table 4: Final results for all octet members. The first two columns give the baryon octet sigma terms 
(Eq. JTJ), the third one the strangeness content of the octet member (Eq. (f3a|l ). and the last two the 
dimensionless sigma terms (Eqs. (|3bjl . fl3c])). The first error is statistical, the second one systematic. 



perform our analysis only with < 550 MeV the systematic uncertainty of a n N would 



have dropped by almost 50%. (see Fig. 3(a) ). 



6 Results 

Following the method described in the previous section, we arrive at the results quoted in 
Table H 

For the nucleon sigma term our result agrees with the "canonical" determination coming 
from 7r — N scattering data [2] of a n N = 45 ± 8 MeV, but with a larger uncertainty, and 
disfavors the larger value ovjv = 64 ± 7 MeV from [3]. It also agrees with recent estimates 
that use lattice data p^5 | [261431] [39] though these assume a particular functional form for 
the light quark-mass dependence while we consider a whole range of forms. Another study 
considers a variety of chiral forms, but with pion masses only down to 300 MeV [32], yielding 
values in the range 42 to 84 MeV, with a statistical accuracy of about 15 MeV. Our value is 
less consistent with the determinations of [331 134] that are close to the higher value (of [3]) 
for the sigma term. 

We believe that the agreement within error of our results obtained from extrapolations using 
very different functional forms enhances the credibility of our final results. 

In Table Owe give the contribution of each source of systematic uncertainty to the final error. 
Although we show the results for a n w, the conclusion that the chiral extrapolation dominates 
the systematic uncertainty of our computation is generally true. 



Source of systematic error 


error on a^N [MeV] 


Chiral Extrapolation: 




- Pion mass range 


9.0 


- Functional form 


5.5 


Continuum extrapolation 


1.9 



Table 5: Contribution to the total uncertainty of the different sources of systematic error. 



To reduce this source of systematic uncertainty we would need to add ensembles at lower 
pion masses. From the value of the uncertainties of a^N at M n = 200 MeV, we can roughly 
guess the gain that would be obtained by having equally precise data at the physical point. 
Such data would imply a reduction of around 50% in the systematic error. On the other 
hand the situation for both (JgsN and yjv is less clear. Our data set has a much smaller range 
of strange quark masses than light quark masses. If we add the fact that the contribution of 
the strange quark to the nucleon mass seems to be small "per se" , it is natural that the final 
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results for both a S sN and un to show large statistical and systematic uncertainties. Only 
by adding simulation points for a wider range of strange quark masses one can increase the 
chances of getting a precise determination of these quantities without making an additional 
assumption which relates m s to m u d dependence as in SU(3) baryon xPT. 

Regarding the other octet members, the values of the sigma terms agrees with other deter- 
minations based on lattice data and a low-energy effective field theory approach [301 131] . 

6.1 %PT consistency check 

Our objective in this paper is to determine the sigma-terms and strangeness contents of the 
octet baryons at the physical point. To this end, we generated ensembles close to the physical 
value of the strange quark mass and as close as possible to the physical light quark mass, 
with a large enough lever arm in the light sector to perform a credible extrapolation to the 
physical point. 

Such a set of ensembles, however, is not ideally suited to determine the SU(3) LECs, 
which are defined in the chiral limit. That would require simulations with a strange quark 
mass approaching the chiral limit. This means that the LECs extracted from our fits have 
uncertainties which we cannot properly control. Nevertheless, it is worth noting that the 
values obtained are roughly consistent with general expectations from phenomenology |21j 
and other lattice studies [30j[31] (see Table [6]) . This further enhances the credibility of the 
results obtained at the physical point. 



M [GeV] 


b [GeV- 1 ] 


b D [GeV- 1 ] 


b F [GeV" 1 ] 


9A 




0.75(15) 


-0.71(24) 


0.103(60) 


-0.359(72) 


0.92(13) 


0.402(90) 



Table 6: Values of the LECs obtained from our xPT fits. Total errors (statistics and systematics included) 
are computed in the same way as for the sigma terms. For the reasons noted in the text, we do not consider 
these errors to be reliable. 



7 Conclusions 

We computed the sigma terms of light octet baryons by studying their mass dependence 
on quark masses through the Feynman-Hellman theorem. We also obtained the strangeness 
content and dimensionless sigma terms for all octet members. We estimated the systematic 
uncertainties associated with the choice of parametrization for this mass dependence by 
considering two very different expansion schemes. We considered regular expansions of octet 
masses, and a singular three-flavor, low-energy effective field theory approach. Within each 
of the approaches we estimated the impact of higher order contributions. Moreover we varied 
the range of fitted quark masses, cutoff parametrizations and excited-state contributions. 

The most interesting quantities are the ones associated with the nucleon. Our final value 
for a nN = 39(4) (±\ s ) favors the "low scenario" with a nN = 45 ± 8 MeV of [H[5], but the 
size of our uncertainties does not allow us to exclude the higher value of [3]. Regarding 
the strangeness content we obtain a value with large statistical and systematic uncertainties 
2/at = 0.20(7)(l 1 y). The situation is a bit better (thanks to correlations between sigma 
terms and the extrapolated value of the mass) for the case of the dimensionless sigma terms 
fudN = 0.021(3)(l4 1 ) and f§ s N — 0.072(29) (itf q ) , that are the quantities of interest for direct 
DM searches. 
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This work shows the difficulties associated with obtaining precise values for these form factors 
in a model independent way. Even with pion masses as low as 190 MeV the mass extrapolation 
dominates the systematic uncertainties of our final results. This explicitly shows the risks 
of using lattice data and only one approach to perform the mass extrapolation. A proper 
estimation of systematic uncertainties should include the uncertainty associated with the 
model that is used, and, at least in our case, this is one of the main sources of systematic 
uncertainty. 

For a n N the only way to increase the precision of our computation without assuming a 
model would be to generate data at lower quark masses. Replacing the extrapolation by an 
interpolation drastically reduces the uncertainty in the final result, as we have explicitly seen 
by considering the values of the sigma terms at M n = 200 MeV. 

For the case of the strangeness content of the nucleon and (JssN, the situation is even more 
delicate. First, Nf = 2 + 1 lattice ensembles usually bracket the physical value of m s 
with a small lever arm. This is more than sufficient to determine physical observables, 
but it is not the ideal setup to estimate how physical quantities change with m s . Second, 
the contribution of the strange quark to the nucleon mass seem to be small "per se". In 
principle the first problem can be solved by measuring the nucleon mass for a larger range 
of strange quark masses. The second problem is solved if we restrict our analysis to three 
flavor B^PT. This low-energy effective field theory approach relates the m U( i dependence 
of octet members to its m s dependence, and thus constrains j/jv- The value of the nucleon 
strangeness content obtained from the covariant B%PT analysis alone has a small relative 
error un = 0.276(77) (±H), illustrating the success of the idea. However, since we do not 
have a large enough range of strange quark masses to establish the validity of SU(3) B^PT, 
this result must be considered model dependent. Note also that a direct calculation of the 
corresponding disconnected contribution suggests a very small un value, barely compatible 
with the B^PT result [55] , 
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Appendix 



A Sketch of three-flavor covariant baryon 

In this appendix, we give a brief outline of the derivation of our chiral fit form (|12p . Working 
to 0(p 3 ), we use the standard octet baryon Lagrangians [351136] 



C$ B = Ti(\B^D^B-M G BB + ^B 1 ^{u^B} + ^B 1 ^[u^B\) (14) 
4/B = Tr(b D B{x + ,B} + b F B[x + ,B} + ...)+boTv(BB)Tr(x + )+... (15) 
with the chiral tensors 



Dfj,B = d^B + p^B] (16) 

r M = ^(d fl + ...) u + u (d II + ...)^) (17) 

u„ = i(ut(^ + ...)u-u(d M + ...)u t ) (18) 
X± = ± u X j u (19) 

X = 2J3 (M + ...), (20) 

where M denotes the (diagonal) quark-mass matrix. To the order we are working, there is no 
contribution from the 0(p 3 ) Lagrangian, C^ B , to the quark-mass dependence of the octet 
baryon masses. The coupling constants D,F of Eq. (|14|1 are related to the actual couplings 
of the nucleon in the chiral limit through D + F = lim u ^ s _».o 9 A arid F/D = lim u ^ s _>.o £> 
whereas Mq corresponds to the chiral limit value of the mass of the baryon octet. The field 
u(x) = y/U(x) describes a nonlinear matrix representation of the (quasi-) Goldstone boson 
fields and B is the matrix- valued interpolating field for a spin 1/2 octet baryon. The next- 
to- leading order low-energy constants &£>, b F , bo govern the leading quark-mass contributions 
to the mass of a spin 1/2 octet baryon, with Ad denoting the quark-mass matrix for three 
light flavors u, d and s. Furthermore, we note that the parameter Bq is a measure of the size 
of the chiral condensate (0|qq|0) (in the chiral limit) |37j 

R . (0|gg|0)m g -X) 

b = ^ ( 21 > 

where the low-energy constant Fq is identified with the value of the octet Goldstone boson 
decay constant in the chiral limit. Finally, for the calculation at hand one needs to know the 
chiral meson Lagrangian up to 0(p 2 ) [37] 

6$ = ^fTT{{d, + ...)U\d^ + ...)U + x ] U + X U ] }. (22) 

Generalizing the SU(2) calculation of [27J along the lines of [31], we obtain the leading one- 
loop result for the mass of a nucleon to 0(p 3 ) 



M N = M - A(b + (b D + b F )/2)Ml - 2(b + b D - b F )M ] 



_ V 9«nMI [ L Ml M a M a M* a] 4 (23) 

^ 4**F* I V 1 4M 2 arCC ° S 2M + 4M L ° g M$ J + 0{p } ' 

To the chiral order we are working here, the meson- loop couplings g-^N, 9kn, 9r]N can be 
expressed in terms of the two SU(3) parameters F and D, as given in Table [7J For generality, 
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9aB 


N 


A 


s 




TT 

K 
V 


\(D+Ff 
l(5D 2 -6DF+9F 2 ) 

UD-m 2 


D 2 

\{D 2 +%F 2 ) 
ID 2 

3 


i(£> 2 +6F 2 ) 
{D 2 + F 2 ) 

Id 2 

3 


f(F>-F)* 

±(5F> 2 +6L>f+9F 2 ) 

^ (j D+3F) 2 



Table 7: Summary of meson-loop couplings g Q s in terms of the low-energy constants D,F. 



we have expressed the Goldstone boson decay constant Fn entering the loop contributions 
via three different symbols F a to account for the individual contribution of the pion, kaon 
and eta-cloud of the nucleon. 

For a fit to lattice QCD data, it is convenient to slightly rewrite (|23p without changing the 
expression at the order at which we are working: 

1. We express the GMOR-masses for Mk, M^ in (|23p as a function of the GMOR-mass of 
the pion and of the mass-parameter M Ss introduced in ^j: 

Ml -> l -{Ml + M 2 ss ) (24) 
M 2 -> l -{Ml + 2M 2 ss ) (25) 

which amounts to a change of variables. 

2. We identify the GMOR-mass M w of (j4]) with the corresponding lowest-lying 0~-boson 
mass M n in each simulation. Possible deviations in the quark-mass dependence from 
the linear GMOR-behavior predicted in @ can only affect 0(p 4 ) corrections f|23j> . 

3. The Goldstone boson decay constants F a in (|23p are identified with the decay con- 
stant Ff hys . Other assignments compatible with SU(3) symmetry might be chosen as 
discussed in Sec. 13.21 

4. With the same reasoning as in the previous item, at chiral order 0(p 3 ) we can identify 
the chiral limit couplings g-wNidKNig-qN with the physical couplings g^N, 9kn, 9r)N, 
resulting in Table [2j The unknown ratio £ of couplings is either treated as an external 
input to our fits or is left as a free fit parameter. 

With a very similar derivation for the case of the octet hyperons, we finally arrive at our fit 
function (fl^D . 
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